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Abstract 

We investigate the large N limit of the expectation value W{X) of a BPS Wilson loop in ABJM theory, 
using an integral expression of the partition function obtained recently by Kapustin et al. Certain saddle- 
point equations provide the correct perturbative expansion of W{X). The large A behavior of W{X) is 
also obtained from the saddle-point equations. The result is consistent with AdS/CFT correspondence. 
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1 Introduction 



Our knowledge on M-theory has become deeper since the discovery of the worldvolume theory on M2- 
branes [1] [2] [3] [4] [5] [6] . The theory turned out to be a superconformal Chern-Simons theory coupled to 
matterqj. The relevance of the Chern-Simons theory in the context of M-theory was already anticipated 
ini. 

It is well known that a gauge theory may be simplified by considering the planar limit. A set of 
M2-branes in C^'/Zfc background is described by ABJM theory |6] with gauge group \J{N)x\J{N) and 
the Chern-Simons level {k,—k). In the large N limit with the 't Hoof coupling A = N/k kept finite, 
ABJM theory is expected to be dual to Type IIA string theory in AdS4 x CP^ with fluxes . A natural 
expectation is that this correspondence might be directly checked by solving ABJM theory in the large 
limit while keeping A finite. Clearly, this must be a difficult problem. Instead, it may be easier to calculate 
a particular observable in ABJM theory in the limit, and to compare the result with the corresponding 
observable in Type IIA string theory. 

In AdS5/CFT4 correspondence, it is already known that there exists such an observable for which the 
above-mentioned line of research is possible. In [9|, the expectation value of the half-BPS Wilson loop 
was exactly calculated. In fact, it is the exact result both for finite N and finite gyj^jN. It turned out 
that the large N limit of the exact result on the Wilson loop reproduces the corresponding quantity in 
Type IIA string theory which was proposed in jlOj jl2j and to which strong evidence was given in |13| . 

After the discovery of ABJM theory, a BPS Wilson loop operator was constructed in [T3] [15] [H] , and 
the pcrturbative calculation of the expectation value V7(A) was performed. Quite differently from the 
case of A/" = 4 super Yang-Mills theory, the BPS Wilson loop preserves at most i of supersymmetry. This 
fact makes it difficult to identify what should be the corresponding object to the Wilson loop in Type 
IIA string theory. However, it seems to be still reasonable to expect that the BPS Wilson loop would 
have a dual string worldsheet in AdS4 x CP'^. To confirm this conjecture, it is necessary to determine 
the large A behavior of W{X). 

Recently, a localization technique was applied to the partition function and W{X) of ABJM theory in 
[T7|, and an integral representations of them were obtained. It was also shown that the integral repre- 
sentation of W{X) reproduces the pcrturbative expansion of W{X) obtained in [14][T5][16] by calculating 
Feynman diagrams. 

The integral representation of the partition function of ABJM theory obtained in [17] looks similar 
to the partition function of a matrix model where the angular variables are integrated out. It is natural 
to expect that the techniques developed for solving matrix models may be applicable to ABJM theory. 
In this paper, we will show that both the pcrturbative expansion of W{X) and the large A asymptotic 
behavior of M^(A) are derived from two saddle-point equations which are obtained from the integral 
representation. The pcrturbative expansion of W{X) derived from the saddle-point equations exactly 
coincides with the one in |17j . including the phase factor due to the framing. The large A behavior of 
W{X) turns out to be 

WiX) - (1.1) 

^ BLG theory was rewritten as the Chern-Simons-matter theory in [?]• 
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for a constant c > whose upper bound is obtaincqj. This result is consistent with the conjecture 
claiming that a string worldsheet in the bulk would be dual to the BPS Wilson loop. The large N 
solution for finite A, however, seems quite difficult to obtain, and it is still an open issue. 

This paper is organized as follows. The localization calculation of [17] is briefly reviewed in section [2l 
The saddle-point equations are derived in section[3] In section^ the perturbative expansion of W{X), up 
to 0{X^), is derived from the saddle-point equations. Section [5] shows a recursive algorithm for calculating 
W{X) perturbatively. Based on some observations on the saddle-point equations for finite N obtained 
in section [6l we analyze the large A behavior of W{X) in section [71 Section [8] is devoted to discussion. 
Appendix 1X1 contains the details of the calculation of W{X) up to order A^"'^. 



2 Localization for ABJM theory 



In suitable situations, the localization is a very powerful tool to exactly calculate some quantities of a 
supersymmetric theory. It was applied to the expectation value of the half-BPS Wilson loop in A/" = 4 
super Yang-Mills theory in [3]. The result is remarkably simple, namely, the path- integral for the Wilson 
loop reduces to an integral for the Gaussian matrix model which can be performed easily. Interestingly 
enough, this exact result enables us to find the behavior of the Wilson loop in the large 't Hooft coupling 
limit which coincides with the one expected from AdS/CFT correspondence. See also [HI [HI [lOjllT] for 
applications to other theories. 

Recently, a similar localization technique was applied to supersymmetric Chern-Simons-matter the- 
ories including ABJM theory in |17| . There is a technical simplification when the three-dimensional 
theories are considered. One can use superfields which are obtained by the dimensional reduction from 
A/" = 1 superfields in four-dimensions. In A/" = 2 gauge theories in three-dimensions, there exists a Wilson 
loop [22j which preserves a fraction of the supersymmetry realized on the superfields off-shell. In the 
case of A/" = 4 super Yang-Mills theory, the supersymmetry which is preserved by the half-BPS Wilson 
loop is always realized only on-shell, and therefore, a sophisticated construction of an alternative off-shell 
symmetry was necessary in [9]. 

ABJM theory has the gauge group \J{N)xlJ{N). According to this, there are two kinds of Wilson 
loops, each of which is constructed from either a gauge field or of the two U(A^) factors. In the 
following, we only consider the Wilson loop for A^. A BPS Wilson loop was constructed and investigated 
in [Ml [H] [HI . For the fundamental representation, the explicit form is 



Wn(C) - ^Ti-nPcxp / ds\ix''A^{x) + \i\MA^ X^{x)Xb{x)\ . (2.1) 
^ Jc ^ -' 

Here X"^ are scalar fields in the bi-fundamental representation of U(A'') xU(iV), and in 4 of SU(4) R- 

^ In some earlier versions of the paper, we claimed that the value 1 17.471 of c would be the exact value. After the 
appearance of 1251 . the arguments were re-examined and then it turned out that this claim was too strong which could not 
be justified by the analysis done in this paper. As explained in section [7] the best we can find is the upper bound on c 
which is compatible with |25| . 
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symmetry. Ma^ is a constant matrix. It turns out that a fraction of supcrsymmetry is preserved iff 



Ma' 



diag(+l, +1,-1,-1) 



(2.2) 



up to an R-symmetry transformation. 

The locahzation performed in [T7] goes roughly as follows. ABJM theory is defined on S^, and the 
contour C of the Wilson loop is placed on the equator of S^. Let Q be a fermionic transformation induced 
by a suitably chosen supercharge which is preserved by the Wilson loop. By definition, is a bosonic 
transformation of the theory. Let V = ip^Qip where ip is the collective notation for the fermions in the 
theory. If V satisfies / d^xQ^V = 0, then the path-integral can be modified to 

^-SABJM-t J d'xQV 



^ABJM 



without changing the value of the partition function. The same is true for the Wilson loop: 



7-1 

'ABJM 



-SABJM-tJ d^xQVT^^^^y 



(2.3) 



(2.4) 



Those quantities are independent of t. If t is taken to be large, then, since the bosonic part of QV is 
\Qip\'^ which is positive definite, the path-integral is localized to a set of field configurations for which 
the bosonic part of QV vanishes. The allowed configurations turn out to be coordinate independent, 
and therefore, the path-integral reduces to a finite-dimensional matrix integral. There could also be the 
one-loop contribution and the non-perturbative contributions to the integrand of the matrix integral. In 
[T7] . the former was calculated explicitly, and the latter was shown to be absent. 
As a result, the partition function of ABJM theory becomes 



Za 



B.JM 



) lli<£ 



,=i Y{^J cosh [7r( 

and the expectation value of the Wilson loop (|2.ip becomes 



(2.5) 



(Wn(C)) 



7-1 

^ABJM 



N 



1 



N 



N ^ 

i=l 



n 



sinh^f 



)] sinli^[7r(0i - 4>j)] 



(2.6) 



HijCOSh [TT{(j), - (jjj)] 

The explicit expression of (Wn(C)) was obtained perturbatively in [T7] using (|2.6p which exactly repro- 
duces the expression obtained in |14||15]p^ by calculating Feynman diagrams. Note that only planar 
diagrams were taken into account in [14| |15j [16] , while (|2.6|) provides the expression exact in both k and 
N. 



3 Saddle-point equations 



In the previous section, we briefly reviewed the derivation of an expression (|2.5p of the partition function 
of ABJM theory in terms of a finite-dimensional integral. This is an exact formula with respect to N 
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and k. If one is only interested in the large N limit, then the necessary information should be contained 
in the saddle-points of the integral. Indeed, it is sufficient to obtain the large N results when one would 
like to compare some results on ABJM theory with the corresponding results on the classical gravity in 
AdS4 X CP'^. This saddle-point approximation would be a more efficient way to extract information on 
observables in ABJM theory than the exact integral expression, at the cost of losing all 1/N corrections. 
Those corrections would possibly be discussed starting with the saddle-point results, but it is beyond the 
scope of this paper. 



The saddle-points of the integral (|2.5p are obtained as the extrema of the effective action 



N 



i<j 



|logsinh^[7r(0i - (j>j)] + logsinh^[7r(0i - 



Vlogcosh2[7r((/),-<^,)]. (3.1) 



ij 

The saddle-point equations derived from S'off are 

— ikipi — coth[7r(0i — (j>j)] — tanh[7r(t/)i — (3-2) 

ik<j)i = coth[7r(^i — — tanh[7r(0j — cj)j)\. (3-3) 

It is important to notice that the coefficients in the left-hand side of (|3.2pp.3p are purely imaginary. 
Therefore, this set of saddle-point equations cannot have a solution in which all and (ja ^re real, 
although it is expected from the derivation of (|2.5p . In addition, one can easily show that a solution in 
which all (pi and (pi are purely imaginary is not allowed. We assume that the integration contours in (|2.5p 
can be suitably deformed away from the real axis so that a complex solution of the saddle-point equations 
provides the dominant contribution to the integral. As we will show below, the result of the perturbative 
calculation of the Wilson loop expectation value seems to support the validity of this assumption. 

It will turn out in the next section that it is convenient to rescale the variables for the perturbative 
calculation. We define Xi and Xi by 

(pi = —Xi, (pi = —Xi, (3.4) 

TT TT 

where A = is the 't Hooft coupling. In terms of these new variables, the saddle-point equations 
([33) (1231) become 

Vx coth[VX{xi — Xj)] — ^ VXtanh[VX{xi — xj)], (3.5) 
-i- ^ VXcoth[Vx{xi — Xj)] — ^ ^ \/Atanh[-\/A(xi — Xj)]. (3.6) 

Let us consider the large N limit of p.5pp.6p . We would like to discuss the limit: 

iV ^ oo, A: ^ oo, A = fixed. (3.7) 



i 


1 






TT 


TV 


i _ 


1 


Xi — 




TT 


TV 
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In this limit, the distributions of Xi and are expected to become continuous. Let p(x) be the distribution 
function of Xi defined on a curve / in C. and let p{x) be such a function for Xi defined on / in C. As 
mentioned above, the curves /, / do not lie on the real axis, nor on the imaginary axis. The functions 
p{x) and p{x) are required to satisfy 



dx p{x) = 1 = J dxp{x). (3-8) 
In the large N limit, Egs. p.Spp.Op are written in terms of p{x) and p{x) as 

dx' p{x')V\coth[y/\{x - x')] - J^dx' p{x')VXta.nh[VX{x - x')], (3.9) 

dx' p{x')y^coth[VX{x - x')] - / dx' p{x')VXt&nh[VX{i - x')], (3.10) 



I 

—X = 

TT 



I 

— X 
TT 



where x £ I and x £ I arc assumed. 

The expectation value of the Wilson loop ()2.6p is 



W{X) = / dx p{x)e^'^'' . (3.11) 



4 Perturbative calculation of Wilson loop 



In this section, we calculate VF(A) perturbatively in terms of small A. By expanding the right-hand side 
of ()3.9p with respect to A, and truncating terms of order A"+^ or higher in the resulting power series, one 
obtains the following truncated saddle-point equation 



I" J/„ XX L \ 



(4.1) 



where 

{{x~xf)r. ^ ^ dX' p^{x'){x - x')\ (4.2) 

{{x-5:f)n = dx' p,,[x'){x - x'f (4.3) 

for a positive integer /c, and the coefficients Cfc, are defined as 

^ oo 

cotha: = - + ^Cfca;^''"\ (4.4) 

^ k=l 

oo 

tanhx = ^tfex^'^-^ (4.5) 

k=l 

The function pnix) is related to p{x) as 

Pnix) = p{x) mod A"+^ (4.6) 
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Pn{x) satisfies the eorresponding equation similar to (j4.1[) . Note that In, In may change by varying n. 

Remarkably, for the determination of Pn{x), it is sufficient to evaluate {{x — x)''')n and ((a; — 5;)'^)„ only 
up to order A"^^ since they are always multiplied by A in (j4.ip . Therefore, assuming that pn-iix) and 
Pn-i{x) are already known, the sum in (|4.ip can be evaluated explicitly. Then, (j4.ip can be rewritten 
in a form of the equation which is familiar in the context of the one-matrix model. (See (jS.Sp below.) 
Namely, the sum in the right-hand side, moved to the left-hand side, is regarded as a contribution from a 
matrix model action which is a polynomial of a matrix of degree 2n — 1. This saddle-point equation can 
be solved by the well-known techniques, and therefore one can determine Pn{x). In this way, p{x) can be 
determined perturbatively in terms of small A. p{x) is determined similarly. 

W{\) is calculated via the formula (13. lip . The perturbative expression of W{X) up to order A" is 

W^(A)==l + ^^^^(i^-)+0(A"+i), (4.7) 

k=l 

obtained by expanding the integrand. Here {x^) are the averages calculated in terms of p{x). It will be 
shown later that {x^) with k odd vanishes to all orders in A, implying that W{\) is indeed a series of A. 
At this order of A, {x^) can be replaced with {x^)n-i since these are always multiplied by A. Therefore, 
knowing pn~i{x), W{X) is determined up to order A". 

In the following, we will show explicit calculations for lower orders of A. 



4.1 0(AO) 

At this order, ()4.ip becomes 



L.^Ux'P^. (4.8) 

TT J, X~ X' 



Note that po{x) decouples from po{x). If the coefficient of the left-hand side is real and positive, then 
this equation can be solved easily. We will use the following trick. We first solve 

'-x^Ux'^, (4.9) 

Jlo X 

assuming n being real and positive, and then set n = iri. Let us define the resolvent 

Rn{x):= [ dx'^^ {xiln) (4.10) 

X — X 



for each order of A. Ro{x) has the form 



Ro{x) = —X y/x"^ — 2k 

K K 
1 K 

X 2x^ 2x^ 



(4.11) 



It is important to note that (x'^)o = for odd k. 

The function po{x) satisfies (|4.9p with k = —iri. Therefore, all the properties of Ro{x), defined 
similarly to Ro{x), are derived easily from those of Rq(x). 
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The averages (x'"")„ are determined by Rn{x) as 



^ oo 

= T + 



(4.12) 



k=l 



Using the expansion (|4.1ip . Vl^(A) can be determined up to order A, as mentioned above. We obtain 

W{\) ^l + TiiX + 0{X^). (4.13) 
Indeed, this result agrees with that of [T7]. 



4.2 O(A^) 

The saddle-point equation (|4.ip becomes 



1^= / dx'^ + A 



i(a; - i)o - (.T - .x)o 



Using the results {x)q — {x)q —Q obtained in the previous subsection, (|4.14p can be written as 

, Pi{x') 











6 TT 





dx' 



X — X' 



This is again the equation of the form ()4.9p with 



--2-A-1 



(4.14) 



(4.15) 



(4.16) 



3 TT 

As a result, Ri{x) has the form (|4.1ip . W{X) is determined up to order A^, and the resulting expression 
is 

1 



W{X) = 1 + TTiA + -tt'A' + OiX-"). 

o 



(4.17) 



This coincides exactly with the result obtained in |17j . To see the agreement with the perturbative 
calculation [H] [15| [16] , it is convenient to pull out the phase factor e^^^ due to the framing jTJ. The 
result is 



W{X) = e 



TTlX 



1 + -n'X' + 0{X^) 



(4.18) 



4.3 0{X^ 



At this order, a generic phenomenon occurs in the perturbative calculation. Namely, the saddle-point 
equation (|4.ip we have to solve becomes more complicated than the Gaussian one (|4.9p . 
Eq. (|4.1[) at this order becomes 



I" Jl2 



dx 



, P2{X') 



+ A 



1 



(X — X)l — [X — Xjl 



+ X' 



1 

45 



{{x-xf)Q + -{{x-S:)^)o 



(4.19) 
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Using the result (|4.1ip with (|4.16p . this equation can be written as 



+ aix 



J12 



I P2{x') 



where 



a, = --A, 

ai = + -A+— A^. 

TT 3 15 



(4.20) 

(4.21) 
(4.22) 



This saddle-point equation can be solved by the well-known technique. The resolvent should have the 
form 

(4.23) 



R2{x) = a2X^ +aix- {l32X^ + /3i) V^?^^^. 
The constants {32, Pi and 7 are determined by requiring 

lim xR2{x) = 1. 

The result is 

TT 3 9 

7 = 2..+ — A-^A^ 
The expansion of R2{x) in terms of provides 



= y + yA--A^ 



'~8 r-^ + ^^ ' 



which are used to determine W{X) as follows, 



W(X) = l + iTTX+-TT'^X+-Tr^\^+0(X^) 

3 6 

1 + ^n^X^-in^X^ + OiX^) 
6 2 



We obtained the perturbative expression which agrees exactly with 



[nQ. 



(4.24) 

(4.25) 
(4.26) 
(4.27) 

(4.28) 
(4.29) 
(4.30) 



(4.31) 



* An argument was given in 1161 showing that W{X) should be a series of A^. The argument was based on the use of 
a regularization of the Wilson loop in which the loop lies on a two-plane. This would be realized by introducing a set of 
concentric loops with slightly different radii, but a regularization like this would break supersymmetry. As pointed out in 
1171 . in our calculation, a non-trivial framing would be introduced as a regularization in which the loops cannot lie on a 
two-plane for preserving supersymmetry. As a result, the argument in 1161 cannot apply to our calculation, allowing the 
appearance of the A'' term in (14.311 1. 
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5 A recursive algorithm 



The calculation of M^(A) shown above can be done systematically. In fact, there exists a recursive 
algorithm to determine the resolvent i?„(a;), i?„(x) from i?„_i(a;). In this section, we show this 

algorithm. We start with the suitable choice of an ansatz for the resolvents at each order of A. 



5.1 Ansatz for the resolvent 

We claim that, for every positive integer n, the resolvent Rn{x) and i?„(x) should have the form 



Rn{x) 


n 






A.— 1 


\k=l 


Rn{x) 


11 

k=l 


\k=l 



^x^ - 7, (5.1) 



fex2^-2j ^/^^^ (5.2) 

where a^, Pk etc. arc polynomials of A of degree n. The orders, i.e. the minimum power of A with 
non-zero coefficient, of and /3fe with fc > 1 arc fc, while those of q;i,/3i and c arc 0. ctk etc. have the 
similar properties. 

We have shown that Ri{x) has the above form, and it is easy to show that Ri{x) can be obtained 
from Ri (x) by 

ak = al, Pk = Pk, 1 = 1*- (5.3) 
Suppose that i?„_i(a;), Rn-i{x) have the claimed form. These forms of the resolvents imply 

{x'')„-i = 0, {t)n-l = 0. (5.4) 

for odd k. 

Let us consider the saddle-point equations at order A". (|4.ip can be written as 



(5.5) 



where 

n 

I 



Y^akx"^-'^--x-J2c,X^{ix-xr-'} +J2t,X'^{ix-ir-'} . (5.6) 



fc=l k=l k=l 



The right-hand side of (|5.6I) is indeed an odd polynomial of x of degree 2n — 1 since a term with an 
even power of x is multiplied by {^x^)n-i or (i )n-i with k odd which vanish as shown in (|5.4p . As a 
result, the distribution function Pn{x) is also symmetric at order A". For this saddle-point equation, the 
ansatz (|5.1[) is the suitable choicqj. By induction, the resolvent Rn{x) has the form (|5.ip to all orders in 



This is the suitable choice as long as we restrict ourselves to a one-cut solution. Although the right-hand side of l|5.6|l 
is a polynomial of a high degree in general, those would not change the structure of the cut drastically since each monomial 
of a high degree is always multiplied by a high power of A, and therefore the effect should not be relevant. 
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A. The ansatz (j5.2p for is also deduced similarly. As a corollary, the eigenvalue distributions are 

symmetric to all orders in A, implying {x^) = (x ) = for odd k. 
The coefficients ak can be written explicitly as 

- -4,1- (5.7) 

TT 

This can be calculated using Rn-i{x) and 

Since the saddle-point equations for pn{x) and pn{x) decouple, they can be solved separately. In the 
following, we focus on the solution for Rn{x). Rn(x) can be determined similarly. 



7 1^ 



21 - 1 
2k -1 



/ ri-l \ / n-1 



5.2 Determination of Rn{x) 



The resolvent Rn{x) is assumed to be of the form (|5.1|) in which au are given as (|5.7|) . Then, /3fc and 7 
are determined by the requirement 

lim xRn(:x) = 1. (5.8) 

a."; -^00 

This is equivalent to the fohowing requirement 



x^ 



lim Rn{x) = 1, (5.9) 



x— f 00 



which turns out to be more convenient. 

The Taylor expansion of (x^ — 7)^^i?„(a;) is 

+a„x2""2 -t-pi7a„x2"-4 _^p272a„a;2"~'^ • • • +p„-i7"^^Q!„ +p„7"a„a;-2 

+a„_ix^""'* +pi7a„_ia;^""'' • • • +p,i_27""^an-i +p„_i7""^a„_ia;"^ 



-ai +pijaix ^ 



where Pn are defined as 



The condition (|5.9p determines /3fe in terms of ak and 7 as 



Let 7 be 



(5.10) 



{l-x)-^=f2pkx'', Pk^^^^^. (5.11) 
^-^ Jn ft! 



/3fe = ^p,_fcy-'=a,. (5.12) 



fc=0 
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7fe with k < n arc obtained from Rn-i{x), and therefore, they are supposed to be known. Since the 
order of with fc > 1 is larger than 1, (|5.12p determines /3j. up to order A" without knowing 7„. The 
remaining unknown constant 7„ is then determined by requiring 

n 

= 1. (5.14) 

This is equivalent to a linear equation for 7„ which can be solved easily. 

This completes the determination of R„{x). By induction, it can be shown that Rn{x) is obtained 
from Rn{x) by ((O)) . 

The expansion of i?„(x) in terms of provides {x^)n via ()4.12|) . The explicit form of them is 



-^g.+fcf+^A, (5.15) 



i=l 

where 

(l-x)^=f:,..^ .. = -^^^- (5.16) 

These quantities are used to determine Uk at the next order by (j5.7p . In this way, Rn{x) can be determined 
recursively. 

Recall that the vev of the Wilson loop W{X) is obtained as 

w{x) = 1 + E + (^-1^) 

We performed the calculation of W{X) up to order A^^, according to the algorithm described above. 
The rcsuh of W{X) if@ 



W{X) = e 



TTlX 



5n^X^ iTT^X^ 29tt^X^ in^X^ Ibln^X'^ itt^A^ 874497r8A8 

1 " 



6 2 120 12 1008 10 362880 

2603in^X^ 344739l7riOAio 1166161i7r"A" , 

+ 0(A^2) 



15120 7983360 3628800 

Further details of the result is summarized in appendix |Al 



(5.18) 



6 Solutions of the saddle-point equations for finite N 

We have shown that the saddle-point equations (|3.2p(|3.3p provide the correct perturbative expansion of 
the Wilson loop W{X). This shows that the saddle-point equations are convenient tools to provide per- 
turbative results of ABJM theory, enabling one to bypass complicated calculations of Feynman diagrams. 

^ In some earlier versions of the paper, the coefficient of A^^ was not correct. After the appearance of 1251 . the calculations 
were re-examined, and then a small mistake was found which affected only the coefficient of A^^. By fixing it, we found 
that the result perfectly coincided with that in 1251 up to the order we calculated. 
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Using the algorithm shown in the previous section, one can systematically calculate W^(A) up to any 
desired order. Other observables in ABJM theory may also be calculable provided that they are given in 
terms of (pi and (pi. 

However, this perturbative calculation, although it makes the calculations much easier, does not 
provide information on the large A behavior of M^(A) which is one of the central issues in AdS/CFT cor- 
respondence. The correspondence between Wilson loops and string worldsheets has been well-established 
in AdS5/CFT4 case [10] [11] [12] , and it would be natural to expect that a similar correspondence holds 
in AdS4/CFT3 case [H] [S] [16] . This conjecture would imply that iy(A) would behave as 

M^(A)-e"^ (6.1) 

with a suitable constant c for large A. This result is expected to hold in the large N limit with A large 
but finite. Therefore, it is natural to expect that the saddle-point equations might be useful to check the 
validity of the conjecture. 

In this section, we solve (|3.2pp.3p for = 3 and arbitrary k numerically. Of course, this investigation 
cannot give any definite results on the large A behavior of W{\). However, we will see that the results 
for A = 3 may provide us a rough idea on the dependence of the distributions of (pi , (pi on the value of 
k, that is, on A. In the next section, we will discuss an ansatz for the distributions for large A based on 
the observation in this section. We will argue that the ansatz seems to be compatible with the behavior 
(IBID . 

One may notice that p.2pp.3p have symmetries. These are invariant under the simultaneous sign 
flips — > {— In addition, (|3.3I) is obtained by taking the complex conjugation of (|3.2p 
accompanied by a replacement (p* pi. It is clear from these symmetries that, provided that {(f>i,<pi} is 
a solution of p.2|)p.3p . then 

{PIP*}, {~4>t,-pn, (6.2) 

are also solutions. We are led to assume that there is a solution of the form {pi,p*}. Under this 
assumption, p.2p and p.3p are equivalent. Therefore, it is enough to solve p.2p only. 

Now let us focus on the case A^ = 3. By the symmetries mentioned above, it is natural to assume 

pi = z, p2^ 0, p3 = -z. (6.3) 

The only non-trivial equation is now 

F{z) := ikz + coth(7rz) + coth(27rz) - tanh(27rRe(z)) - tan(27rlm(z)) - tanh(7rz) = 0. (6.4) 

The graphs for Rei^(2;) = (funl) and Imi^(z) = (fun2) are plotted in Figll]and Figl^J Here we focus 
on the solution which is the closest to z = 0. 

Note that the solution is almost on the line y ^ x for k = 100 where z = x -\- iy. This seems to be 
consistent with the perturbative result obtained in sectional In fact, the solution at order A*^ describes 
a distribution of pi which are on the line y = x. On the other hand, for fc = 1, the value of x increases 
compared with the k = 100 case, but the value of y does not increase as much as x, and therefore the 
solution z is placed apart from the line y = x. 
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Based on this simple observation, it is tempting to speculate that the distribution of 0^ may become 
broadly extended in the x-direction, while in the y-direction the width of the distribution would be of 
order A'' in the large N limit with A large but finite. In the next section, we will argue that the appropriate 
rescaling of the variables would be 

7r0i = A{X)xi + if{xi), (6.5) 

where A(A) is an increasing function of A, and f{x) is a function of order A*^. The new variables Xi are 
assumed to have a distribution function p{x) in the large TV limit which is another function of order A'^. 
These two unknown functions p{x) and /(x) are expected to be determined by the real part and the 
imaginary part of (|3.2|) . The function A(A) will be determined by the consistency of the order of A. 

The assumption (|6.5I) seems to be appropriate by the following reason. Notice that the right-hand 
side of p.2p has the term 

-tanhK0, -0*)], (6.6) 

since the second sum does not exclude j — i- Due to this term, lm(0i) cannot approach ±^ with a finite 
value of Re((/>i) since, if this is the case, the right-hand side diverges while the left-hand side is finite. In 
the large A'^ limit, (f>i are expected to form a continuous curve / in C. The above observation suggests 
that / should be confined in a region {z £ C\ — j < Im(z) < 



7 Toward large A 



In this section, we try to extract information on the distribution of for large A as much as possible. 
With the ansatz (pi — (j>*, the saddle-point equations p.2|)p.3|) reduce to 



- l^'^ = ^ 5Z coth[7r(0, - 0,)] ~^Y1 tanh[^(0, - cj)*)]. (7.1) 

This equation is still complicated. Let us start with the discussion on a simpler equation, and then 
gradually increase the complexity of equations. 



7.1 coth-model 

The first equatior0 we consider is 

where c > is a constant of order A'^. The distribution of (f)i is symmetric with respect to the origin 
(j) = due to the symmetry of the equation. Suppose that the width of the distribution is proportional 

' This equation is actually the saddle-point equation for a matrix model describing Chcrn-Simons theory on [26] . 
This matrix model was further investigated in 1271 . See also 1241 . 



15 



to A(A) which satisfies 

lim A(A) = oo, (7.3) 

A— f oo 

and the proportionahty coefficient is of order A*^. This assumption seems natural wlien (j7.2p is regarded 
as an equation for the balance between an external confining force acting on 0j (left-hand side) and a 
repulsive force among (pi (right-hand side). As A becomes large, the external force which confines <f>i 
around = becomes small while the repulsive force is kept intact. As a result, the distribution of 4>i 
becomes broader as A becomes larger. 
Define rescaled variables 

x^ := A^V.- (7.4) 

By definition, the width of the distribution of Xi is of order A° . Consider the large N limit while keeping 
A large but finite. Then (|7.2p can be written as 



cA 



dx' p{x')coth[A{x-x')l (7.5) 



in terms of the distribution function p{x) of Xi. We denote the support of p{x) as /. As mentioned above, 
/ should be of the form 

/ = [-.T,„, Xrn] C R, X„, = ©(A^). (7.6) 

Ea. (|7.5p can be rewritten as follows, 
where e{;x) is the sign function. 



Notice that the integrand of the second term in the right-hand side of (|7.7|) is localized in a region 
around x with width of order A^^. Therefore, only some local information on p{x) determines the 
dominant part of the integral for large A. Based on this observation, let us assume that there is a 
neighborhood U {x) C / of a; in which p{x) is smooth, and its Taylor expansion at x has coefficients of 
order A". Under this assumption, the integral is estimated as follows. First, consider an integral 



i 



b ^n^-A\x\ r+oo ^n^-A\x\ r-a ^n^-A\x\ r+oo x"e~^ 



dx . ^ ^ — =f dx . ^ ^ / drz . ■ . / dx^——, (7.8) 

smhAa; ./-oo smhAx ./_^ smhAx smhAx 



where a < < 6 and n> 0. The third term in the right-hand side is estimated as 

dx ' . - &"A-ie-2A6 ^ o(A-2e"2Afc)^ (7 9) 

sinh Ax 

The second term is also estimated similarly. Both of them are negligible compared with the first term 
for large A. Therefore, the integral (|7.8p is 



J a 



■^^2.«g-A|x| I A-"-i7i!2-"C(" + l) (?^:odd), 
sinh Ax I 0, (n : even). 



This estimate allows one to obtain 



16 



This term is actually negligible compared with the first term in the right-hand side of (|7.7p . For large A, 
p{x) satisfies 

cA f 

—x^jdx'p{x')e{x~x'). (7.12) 

This implies 

A(A) = A, p{x) ^ \. (7.13) 

The above analysis fails where x is very close to one of the boundaries of / so that the width of U {x) 
is less than or of the same order of A^^, and where p(x) starts changing rapidly. The simplest possible 
solution for p{x) would be such a fmiction which is almost constant except for regions jx — Xm| = 0(A^^) 
and \x + Xra\ ~ 0{Sr^), and which rapidly decreases to zero near the boundaries of /. For this solution, 
Xra is determined to be 

x^ ^ (7.14) 

for large A. 

The results obtained so far turn out to be enough to deduce an approximate solution of (|7.5p which 
may converge to the exact solution in the large A limit. The approximate solution is 



■ tanh[A(x + c^^)] tanh[A(x - c^^)], (|x| < c^^) 
0, {\x\>c-^). 



Pi^) = S „ /, , , (7-15) 



Let f{x) be defined as 



f{x):=4- dx' pix')coth[A{x^x')]. (7.16) 



If X is not equal to ±c ^, then the argument above is valid for a large enough A, and therefore, (j7.5p is 
satisfied. On the other hand, if a; = +c^^; then /(c~^) becomes 

/(c"^) = ^ t dx' tanh[A(x' + c'^)] (7.17) 

which satisfies (|7.5p in the large A limit. The case x = ~c^^ is similar. 

The graphs of f{x) with c = 1, which is expected to be a good approximation to x, are plotted in 
Figl3] for A = 50, 100. The graphs indicate that (|7.15p would be a good approximation in the large A 
limit. 

We have found that the distribution of 0i satisfying (|7.2p has a width 2c" ^ A for large A. This width 
is larger than the one in a model where cothx in the right-hand side is replaced with x~^, that is, in the 
case of the Gaussian matrix model. This is simply because cothx provides an infinite-range repulsive 
force. 

In fact, the above information on p{x) is enough to determine the large A behavior of the "Wilson 
loop" 1^1 (A): 



Wi{X) := / dxp{x)e^'' 



e 



(7.18) 
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Figure 3: Plot of f{x) with A — 50, 100. f{x) approaches to a; as A becomes large. 

A derivation of this estimate is as follows. Consider the estimate of 

W{K)^ r ^ dxp{x)e^-' (7.19) 

for large A. This can be written as 

M^(A)=e^^'"/ -u)e"^". (7.20) 

Since A is large, only the information of p{xm — u) around u = is important. Suppose that p{xm, — u) 
behaves as u"p(u) where p{u) is regular at u = 0. Then, the leading term of the asymptotic expansion 
of W{X) with respect to A~^ is 

/"OO 

1^(A) - 6^=^- / duM"p(0)e-^" 
Jo 

= p(0)r(a + l)A-"-ie^"'". (7.21) 



It should be noted that the coth- model was solved exactly in the large N limit in [IZHH]. The 
resulting distribution function has the behavior which we found in the above analysis. 



7.2 cosech-model 



Next, let us consider 



c 



-77 coth((/), - 0j ) - ^ XI tanh((/), - ) 



-Y • (7.22) 

iV^sinh[2((/), -(/.j)] ^ ' 
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In this equation, the infinite-range repulsive force due to cotha; terms is canceled by tanhx terms. This 
equation looks closer to (j7.ip . 

As in the previous subsection, define rescaled variables 

X, := A- V^ (7.23) 

so that always of order A*'. In the large limit, (|7.22p becomes 

cA f , , 2p(x') 

—X = 4- dx' . , , ; ' —. (7.24) 

A J sinh[2A(a; - x')] 

As in the coth-model, there is an integral in the right-hand side whose integrand is localized around x 
for large A. Therefore, the integral can be estimated similarly. One obtains 



A difference from the coth-model is that there is no 0(A°) term in the right-hand side of (|7.24p . As a 
result, the right-hand side must be small for large A, implying a different functional form of A(A). The 
appropriate choice of A in this case is 

A(A) = A3. (7.26) 

This is a reasonale result. The width of eigenvalues is smaller than that of both the coth-model and the 
Gaussian model, refiecting the fact that the range of the repulsive force among the eigenvalues in the 
cosech-modcl is the shortest. 

The equations (|7.24p and (j7.25p determine p{x) as 

pix)^^ixl-x^) (7.27) 
where Xm is an integration constant. It turns out that a;„i is determined by requiring 

dx^{xl-x')^l. (7.28) 

The analysis of the coth-model suggests that an approximate solution for the large A limit would be 
obtained by modifying p{x) slightly at the boundary of its support such that the modified distribution 
function satisfies (|7.24p at x = ±x,n- A possible modification is 



2c 

P[X) = 2 



2 „ 2 I _£ 



tanh[2A(x - Xm)] tanh[2A(2; + Xm)], (7.29) 



where 

^^TT--y . (7.30) 

^ TT ^ (2n -f 1)2 ^ ^ 

This function p{x) satisfies the correct normalization condition due to (|7.28p in the large A limit. Note 

that the addition of the constant since this is small, does not contradict with the result (|7.27p . The 
constant ^ is chosen such that lim;j_>.oo g{xm) — c.t,„ is satisfied, where g(x) is defined as 

g{x) := A^ r^'^dx' . , y}""^ ^. (7.31) 

J-x^ sinh[2A(x - x')] 
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Figure 4: Plot of g{x) with A — 50, 100. g{x) approaches x point-wise. 

The graphs of g{x) with Xm = 1 and c = 1, arc plotted in Fig|3]for A = 50, 100. By the argument 
using the local nature of the integrand, one can show that g{x) converges to x point-wise in the large A 
limit. 

The large A behavior of the Wilson loop W2{\) in this case is 

W2{\) r dxp{x)e^^'' 

~ e^'"^*. (7.32) 



The FigUl however, suggests that the maximum deviation of g{x) from x, or sup^fz[-i^i]\g{x) — x\, 
might not vanish in the large A limit. Indeed, it can be shown that 

limg(x) = Ar ^^ tanh(t^-y) 
A^-oo TT^ J-oo smhit 

holds, where y = 2A(a;m — x) is assumed to be of order A°. This implies that 

lim sup \g{x) - a;| > 0. (7.34) 

As a result, the convergence of g{x) to x is not uniform. This may not be a serious drawback if one is 
only interested in the large A limit. If one would like to discuss also a 1/A corrections, it might be better 
to start with another approximate solution. 
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7.3 ABJM equation 



The experiences in analyzing the previous models will be helpful to consider the saddle-point equation 
for ABJM theory. Recall that the saddle-point equation we would like to study is 

1 ^ coth[7r(0, tanh[^(<^, - 0*)]. (7.35) 

In the following, we assume the existence of a solution compatible with the following scaling behavior 



(t>i = -[Axi + if{xi)], 

IT 

where € R and f{x) is a real function of order A''. In the large N limit, (|7.35p becomes 

dx' p{x') coth[A(a; - x') + i{f{x) - f{x'))] 
dx' p{x') tanh[A(a; - x') + i{f{x) + fix'))]. 



(7.36) 



A i , , 



(7.37) 



This consists of two independent equations which determine two real functions p{x) and f{x). 

As in the previous subsections, the integrand in the right-hand side is localized around x for large A. 
Explicitly, 

1 1 



coth(a + i5_) = tanha + 
tanh(a-fi6+) = tanha + 



cosh^ a tanli a + i tan 6_ ' 

1 i tan 

cosh^ a 1 + « tanh a tan b+ 



(7.38) 
(7.39) 



where a = A{x — a;') and b± ~ f{x) ± f{x'), and the infinite-range part tanha cancels, as in the cosech- 
model. Therefore, assuming the smoothness of p{x) and f{x) in a neighborhood U{x) C R of a;, only 
local information on those functions determines the integrals. 

First, we separate the complex equation (|7.37p into its real part and imaginary part. It is possible to 
show that the expansion in terms of 6_ provides an asymptotic expansion of the corresponding integral 
with respect to A^^. In addition, the expansion in terms of (tan 6_/ tanha) also provides sub- leading 
terms in A^^, due to the locality of the integrand. Taking into account only the leading contributions to 
the asymptotic expansion, (|7.37p is approximated as 



A 



dx' 



2p{x') 



1 



tanh[A(x - a;')] tan2(2/(a;)) 



dx' p{x') 



sinh[2A(a; - x')] coth^[A(x - x')] 1 + tanh^[A(a; - x')] i&ii^{2f{x)) 
fix) -fix') , 1 tan(2/(x)) 



(7.40) 



sinh^[A(a; - a;')] coth^[A(a; - a;')] 1 + tanh^[A(a; - a:')] tan2(2/(a;)) 



4A-V(x)/(a;). 



(7.41) 



Here we assumed |/(a;)| < j. Note that naively the order of A of the second integral in (|7.40p is A ^, 
but the contribution of this order vanishes. These equations are consistent with the scaling assumption 
I^LM iff 



A(A) = VA. 



(7.42) 
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p{x) and f{x) are determined by the above equations. The result is 

p{x) ~ -^V2c2 -a;2, (7.43) 

where c > is an integration constant. 

Note that the saddle-point equations are consistent with our ansatz for f{x) to be confined within 
[—J, +f ]• In the case j < f{x) < ^, for example, (|7.4ip is replaced with 

■^x ^ A-^p{x)(4:f -2tt). (7.45) 
ttA 

This implies p{x) < 0, strongly suggesting that the solutions (|7.43p(|7.44p do not extend to a region where 
I /(a;) I > J, that is, p(x) is non-zero only if < c. Since p(c) > 0, this implies that p{x) suddenly starts 
decreasing to zero at a point in the interval [0,c], just like the coth-model discussed before. 

Let Xm < c be the position of the right-edge of the support of p{x). Xm is determined by requiring 

dxp{x) = l, (7.46) 



provided that c is given by another input. Since we do not know at present how to determine c, the 
exact value of a;,„ cannot be obtained from the analysis explained so far. The best we can do is to put 
an upper bound on Xm which is derived as follows. The approximate form (j7.43p of p{x) indicates that 
Xm becomes large if c becomes small. Therefore, under the condition Xm < c, the largest value of Xm 
determined from (|7.46p is obtained when c is chosen to be equal to Xm. For this particular choice, (|7.46p 
determines c to be 

(7.47) 



TT-f 2 

This turns out to be the minimum possible value of the integration constant c suitable for our proble: 
This gives the upper bound 

Xm < (7.48) 
The behavior of H^(A) in the large A limit is then expected to be 

|M^(A)| - e^'^^"^. (7.49) 

This exponential behavior and the upper bound (|7.48p on the coefhcient are consistent with the con- 
jectured correspondence between the Wilson loop and a string worldsheet in AdS4 x CP"^ proposed in 
[SI [E] [E] ■ Especially, the appearance of \/A in the exponent strongly suggests the existence of a dual 
string worldsheet in AdS^. 



^ In some earlier versions of the paper, we claimed that the value l|7.47p of c would be equal to Xm- After the appearance 
of I25| . the arguments were re-examined and then it turned out that this claim was too strong which could not be justified 
by the analysis done in this section. As explained above, the best we can find is the upper bound on Xm which is compatible 
with [25] . 
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8 Discussion 



We have investigated the saddle-point equations (|3.2p(|3.3|) . For small A, wc found that the equations 
reproduce the perturbative expansion of the expectation value iy(A) of the Wilson loop (|2.ip . We also 
found an efficient algorithm to determine the perturbative expansion to any desired order of A. For large 
A, we found an approximate solution to the saddle-point equations which is expected to converge to the 
exact solution in the large A limit. Based on the solution, we determined the large A behavior of Ty(A) 
which may be compatible with AdS/CFT correspondence. 

The algorithm found in section[S]is rather simple. This could provide a recursive formula for {x'^), and 
therefore for W{\). It is very interesting to analytically solve this recursive formula to obtain a closed 
formula for W{X). It is already interesting if one can determine all the coefhcients of the perturbative 
expansion using the recursive formula since some non-perturbative information may be extracted from 
them. 

The understanding of the solution to (|3.2p ()3.3|) for large A obtained in this paper is still at the primitive 
level. Especially, the explicit form of p{x) may have ambiguities at finite A which cannot be fixed by the 
argument of locality and the large A limit. For example, in the cosech- model, although the expression 
(j7.29p provides an approximate solution which converges to the exact solution in the large A limit, it is 
not clear whether it is appropriate to use the expression to discuss 1/A corrections. It is quite important 
to obtain the solution for finite A which can interpolate the small A result and the large A result. One 
may expect that the holomorphicity of 



as well as the conditions derived from the saddle-point equations might be enough to determine p{x) 
uniquely. Indeed, TZ{z) turns out to satisfy some equations which are similar to the one discussed in 
[23] . In the case of ABJM theory, the periodicity of Tl{z) in the imaginary direction makes the problem 
complicated. It is very interesting to solve Eqs. p.2|)(|3.3p exactly, for example, in this way. Another way 
to gain some information on the large A solution is to solve (|3.2|)p.3|) numerically for a large N, extending 
the analysis of section [51 

Since the localization of [17j can be applied to a quite general family of supersymmetric Chern-Simons- 
matter theories, one can extend the analysis exhibited in this paper to more general theories. A systematic 
research in this direction may shade some light on the understanding of AdS4/CFT3 correspondence. 
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Appendix 



A Calculation at order A 
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The result is the foUowing: The parameters in the resolvent Riq{x) are 



ai 

0-2 

as 
a4 

OL'J 

as 
ag 

/3i 



95034445994244677 i tt^ 



824385881240071 7r*A^ 1792539152 i tt^ A« 8481442 tt^A^ 



+ 



11693788644037500 155917181920500 
48010852 iTr^A^ 34774 tt^A^ 2288 ^tt^A^ 6047r2A3 



516891375 
8i7rA2 



3648645 



2A 

30405375 31185 2835 945 ^ 15 ^ T ' 

1300145665611332 TT^A^" 344633417596672 ivr^A^ 186341275429 77*^ A^ 



5 \7 



27842353914375 12993098493375 
4084264 TT^A^ 239872 ^tt^A^ 3907 tt^A^ 128 ^ttA^ 



12405393000 
14 A2 



7321664 iTT^ A 
868725 



868725 



93555 



2835 



189 



45 



441275989540096 iTT^A^o 75094926421943 tt^A^ 183862829632 iTr^A^ 64951384 tt^A^ 



4218538471875 
27517696 iTT^A" 482 7r2A5 



1484925542100 
128i7rA4 124 A^ 



7753370625 



6081075 



6081075 275 225 945 

2535162760294933 tt^ A^o 87217425750016 i tt^ A^ 



151498931417 TT^A^ 752263168i tt^ A^ 



21417195318750 
73397228 TT^ A^ 4096 ^ttA^ 



1856156927625 
254 A* 



8683775100 



+ 



127702575 



42567525 10395 
30838712389203968 i tt^ A^" 
417635308715625 
2830336 iTrA^ 292 A^ 



4725 ' 
4653424941601 TT^A^ 
202489846650 



11226873856 iTT^A^ 50380768 tt^ A^ 



1776226725 



34827975 



11609325 13365 ' 
6113412629359571 Tr^A^" 

235589148506250 
353963361107968 iTT^A^o 
70676744551875 
927244437672232 TT^ Ai° 



230520782848 iTT^A^ 43219333 tt^ A^ 65536 ^ttA^ 5657908 A« 



38940355125 
4193522922926 TT^ A^ 



5568470782875 



39760875 
118521856 iTT A® 
1550674125 



467775 
65528 A^ 
18243225' 



638512875 



91995766784 iTrA-^ 



1891642280653125 
183092903936 i7rAi° 

8955143071875 
366185109428 A^° 
^1531329465290625' 
25177793483209721 i tt^ A^^ 



2292899734125 
22998766228 A^ 
38979295480125' 



33862354 A** 
23260111875' 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 
,(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.IO) 



9640092302437393 A^ 3455382824 ivr^A^ 3032042257 tt^A^ 



20581068013506000 6236687276820000 
3237679 iTT^A^ 436883 tt^A^ %Tli-n'^\^ 346 7r2A3 



9823275 935550 
1308690820561018067 tt^ A^o 



2835 945 
744839342044528 iTT^A^ 



68603560045020000 64965492466875 
372482207 TT^ A^ 734336 ^tt^A^ 4618 tt^A^ 172 ^ttA^ 



5746615875 

2 i TT A^ 2 A i 

6874823870267 TT^ A^ 
976924698750 
14 A2 



3831077250 

(A.ll) 
128817718 iTT^A^ 



30405375 



141891750 



467775 



4725 



315 



45 



(A.12) 



24 



17639281196348848 iTT^A^o 88912240611058 tt^A^ 186517708697 iTr^A^ 1542519221 tt^ A^ 



306265893058125 2998407344625 
2134621256 iTT^A^ 2248 tt^A^ 2434 ^ttA^ 124 A^ 



12524675625 



638512875 1575 4725 945 ' 

5091572820200602367 tt^ A^^ 2166511481052202 i tt^ A^ 



4 \8 



64315837542206250 64965492466875 
963077002 A^^ 6964 ^ttA^ 254 A^ 



12902809900867 TT-* A 
976924698750 



638512875 18711 
607486105446138962 i tt^ 
10719306257034375 
50003524 iTrA^ 292 A^ 



4725 ' 
1210297833352637 TT^A^ 
64965492466875 



3 \8 



2651143334458 iTT^ A 
488462349375 



212837625 

(A.13) 

3053312744 iiTT^ A^ 
^ 638512875 

(A.14) 

280467644 TT^ A^ 



212837625 



212837625 13365 ' 
26234503780513031 tt^ 
1191034028559375 
5657908 A^ 



60442600532504 i A^ 23611974296 tt^ A^ 355768 ^ttA^ 



11464498670625 



23260111875 



2606175 



638512875 ' 

147061705524924976 i tt^ A^o 



139315307981924 TT^A^ 1743965486 i ttA^ 65528 A^ 



32157918771103125 



19489647740062f; 



23260111875 



18243225' 



201731372996882 tt^A^o 
428772250281375 

30942701463628 ittAI" 
1531329465290625 
366185109428 Ai° 



1580440276 iTT A-* 

39978764595 
22998766228 A^ 
38979295480125' 



33862354 A^ 
23260111875' 



1531329465290625 
161894160651233 i7r"Aio 

24198786612000 ^ ~ 
5053207 ZTT^A^ 13922 tt^A^ 



535693251913 TT^^A^J 10384955081 i tt'-' A» 60232 tt^A^ 



3308104800 



109994484600 

2977i7r5A'^ 296 7r''A3 SAin^ An^ X 



3153150 
can be expanded as 



10395 



3150 



Rioix) 



1 

X 



k=l 



315 



2fe\ -2k- 
10^ 



45 



25025 
2 in. 



(A.22) 



(A.15) 

(A.16) 
(A.17) 
(A.18) 
(A.19) 
(A.20) 

(A.21) 



22485647 iTT^A^o 227813 tt^OA^ 3307813 i tt^ A*^ 863 7r8A^ 961i7r^A6 177r*^A5 



49397040 
283i7r5A4 TT^X^ 



453600 
in^X^ TT^X 



2700 6 
12000329219 ttI^AI" 



9 3 
276624839 i tt" A^ 



13608000 3024 

i TT 

T' 

228148609 Tfio A** 



6615 



90 



68188 iTT^A^ 7799521 TT^A^ 



4322241000 
158i7r^A^ 1469 TT^ A'' 



225 2700 
5323721653441 iTfi^Ai" 



130977000 

8 i TT^ A^ 

15 ^ 60" 
413232004 TT^^A^ 



158760000 
297r4A2 2iTT^X 



59535 



9525600 



449513064000 49116375 
2267093 i TT^ A^ 74381 tt^A^ 4811«7r^A'^ 



TT 

37102261243 i7r"A^ 
6286896000 ^ " 
IOItt'^A^ AUtt^X^ 



3448213 tt^OA^ 



793800 
5^4 A 



5 i TT'^ 



714420 30240 
5545927036934549 tt" A^o 

121368527280000 8756748000 

224690939 TT^o 19561 iTr^A^ 51349 tt^ A"* 



2520 63 30 

275906240221 i tt^^ A^ 44102510327 tt^^ A^ 



2020788000 
1418i7r^A3 203 7r6A2 



13002184 i7r"A^ 
841995 



20412000 



2430 



8640 



315 



60 



(A.23) 



(A.24) 



(A.25) 



(A.26) 



25 



10\ 
AO 



12\ 



14\ 



16\ 

no 



;io 



20 \ 



3376146566293313 iTri^Aio 2338727068919 tt" 14226515947H tt" A*^ 4681993579 tt^^A^ 

20228087880000 20803910400 1868106240 89812800 

33573821 i7r"Af5 16586929 tt^o A^ 523711 ^tt^A^ 24103 tt^A^ Giin'^ idn^ X 



935550 665280 30240 2016 8 8 

+ (A.27) 
24761546379007051 ttI^A^" 3010228043441 iTr^^A^ 1081813416673 tt" A*^ 186678329 iTr^^A^ 



10 ^ .oo...n...onnnn + 



42247941120000 7801466400 4245696000 1105650 

7915650007 Tfi^A^ 307303 i7r"A^ 9712351 tt^o A'* 5753 ^tt^A^ 1133 tt^A^ 33i7r^A 

70761600 ^ 4158 201600 189 ^ 64 ^ 4 

33 TT^ 



(A.28) 

9588084588732887 iTT^Ai" 1547644701254461 tt^^ A^ 5920625808253 iTri^A^ 23723634703 tt^^A^ 



4809475440000 1207084032000 7185024000 44906400 

818216699 iTfi^A^ 1807931957 tt^^A^ 16755713 ^tt^A^ 324181 tt^^A^ 7007 ^tt^A^ 



2432430 8553600 129600 4320 180 

1001^ _42W 

64 128 ' ^ ' 

7920651425068514707 tt^^ X^° 20720319106004807 i tt^^ A^ 1461510279337891 tt^^ A^ 



1200445069824000 5001854457600 564583219200 

223063873 iTfi^A^ 1336809970387 tt^^A^^ 41814259i tt^^ A^ 34944793 tt^^ A" 19487 iTr^A^ 



138996 1362160800 71280 103680 108 

40469 tt^OA^ 715i7r9A 715 tt^ 



480 24 128 

87296746879362569 i tt^^ X^° 667329042674322019 tt^* A^ 5221907328724763 i tt^^ A® 



(A.30) 



4095463680000 51111386726400 658680422400 

92034765828557 Tfi^A^ 562945732229 iTr^^A*^ 60332915 tt" A^ 



19372953600 201801600 38016 

3952891 i 7ri3 A^ 4899791 tt^^ A^ 115583 iTr^A^ 7293 7riOA 2431 ^tt^ 



4608 11520 640 128 256 ' 

6470504221105366364567 tt^" A^o 611939373233888527 i tt^^ X^ 337680661534598713 tt^^ A* 



(A.31) 



96035605585920000 15205637551104 14227497123840 

1248183306883 iTT^^A^ 60061822300489 ttI^A^ 558837481 i tt^^ A^ 



90810720 7749181440 133056 

517647229 ttI^A-* 994517 itt^A^ 1767779 tt^^ A^ 20995 ^tt" A 4199 tt^o 



241920 1008 4608 192 256 

The Wilson loop is therefore 

116616H7r"A" 3447391 tt^o A^o 2603 ^tt^A^ 87449 tt^A^ ^tt^A^ WItt^A^ 



(A.32) 
(A.33) 



W{X) 



3628800 7983360 15120 362880 10 1008 

iTT^X'' 297r4A4 in^X^ Stt^A^ , 
-—2 ^0 ^ + ^ + (^-^^^ 
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